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ABSTRACT 


This  report  describes  the  theoretical  performance  of  an  adaptive 
antenna  array  in  a  multipath  signal  environment.  The  array  response 
to  a  coded  communication  signal  comprised  of  a  direct  (line-of-sight) 
component  and  a  specular  multipath  component  is  determined  for  the  case 
in  which  differential  doppler  frequency  shift  between  components  is 
negligible  compared  to  the  reciprocal  time  constants  of  the  array 
processing  loops. 

Performance  analyses  are  presented  which  determine  the  ratio  of 
powers  in  the  direct  and  multipath  signal  components  at  the  array  output. 
The  fraction  of  total  power  at  the  array  output  which  is  usable  or 
beneficial  for  post-array  data  demodulation  is  also  calculated. 

Computer  evaluations  are  presented  which  illustrate  the  dependence  of 
these  ratios  on  received  signal  amplitudes,  carrier  phases,  angle-of- 
arrival  separation,  and  relative  synchronization  with  a  coded  reference 
signal  which  is  locally  generated  at  the  receiving  array. 
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I.  INTRODUCTION 


Adaptive  antenna  arrays  have  been  under  study  for  several  years. 
Widrow,  et.al.[l]  suggested  the  basic  feedback  algorithm  and  did  com¬ 
puter  simulations  of  the  behavior  of  such  an  array.  Shor[2]  and 
Applebaum[3]  have  also  studied  closely  related  ideas.  Baird  and 
Zahm[4],  and  Baird,  et.al.[5]  have  investigated  the  relationships 
between  different  performance  criteria  for  adaptive  arrays.  Early 
experimental  measurements  of  array  performance  in  a  communications 
environment  are  described  in  Riegler  and  Compton[6].  Compton[7]  also 
discusses  the  results  of  extensive  measurements  with  a  four-element 
array  in  the  200-400  MHz  band  where  the  array  elements  are  mounted  on 
an  irregular  surface.  The  use  of  adaptive  arrays  for  clutter  rejection 
in  radar  has  been  studied  by  Brennan  and  Reed[8].  Berni[9]  discusses 
techniques  which  employ  an  adaptive  array  for  angle-of-arri val 
estimation. 

The  application  of  adaptive  arrays  to  interference  rejection  in 
coded  communication  systems  has  recently  been  under  investigation,  both 
theoretically  and  experimentally.  The  behavior  of  an  array  with  coded 
signals  and  a  "bootstrap"  reference  signal  generation  circuit  has  been 
studied  with  computer  simulations  by  Reinhard[10,ll] .  Two  versions  of 
a  coded  adaptive  array  have  also  been  implemented,  and  the  experimental 
results  are  described  by  Huff  and  Reinhard[12,13],  and  by  Schwegman 
and  Compton [14]. 

The  purpose  of  this  report  is  to  investigate  the  response  of  an 
adaptive  array  to  a  coded  communi cation  signal  when  a  specular  multi - 
path  component  is  also  present  in  the  received  signal.  Of  particular 
interest  is  the  case  where  the  direct  signal  and  the  multipath  signal 
arrive  with  only  a  small  timing  difference  (relative  to  the  reciprocal 
of  the  signal  modulation  bandwidth),  because  this  case  appears  to  be 
potentially  the  most  descriptive  to  array  performance.  In  this  report, 
we  study  the  array  response  to  such  a  correlated  interfering  signal 
from  a  theoretical  standpoint. 

In  Section  II  analytical  models  of  the  adaptive  array  processor 
and  the  input  signal  structure  are  defined.  Certain  restrictions  in 
the  models  which  are  necessary  in  order  to  obtain  analytical  solutions 
are  discussed.  In  Section  III  the  performance  resulting  from  these 
models  is  analyzed  for  the  case  in  which  the  array  response  has  reached 
steady-state  conditions.  Ain  expression  for  the  ratio  of  powers  in  the 
signal  components  at  the  array  output  is  derived  and  evaluated 
numerically  in  Section  III-C.  In  section  III-D  the  power  contained  in 
that  portion  of  the  composite  array  output  signal  which  is  correlated 
with  the  array  reference  signal  is  compared  to  the  power  present  in  the 
undesired  output  components.  T’°  numerical  results  presented  yield 
considerable  Insight  into  the  array  behavior. 
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II.  ANALYSIS  MODEL 


A.  Array  Processor  Model 

The  general  configuration  of  an  m-element  adaptive  array  is 
shown  in  Fig.  1.  The  receiving  electronics  for  each  array  element 
which  normally  precede  the  adaptive  processor  are  not  shown  for  sim¬ 
plicity.  These  electronics  include  an  rf  pre-amplifier,  a  mixer  for 
down con vers ion  to  an  intermediate  frequency  (IF),  and  an  IF  bandpass 
amplifier  which  rejects  the  upper  sideband  product  of  the  mixer  output.* 
The  array  processing  consists  of  separating  the  composite  IF  signal 
plus  thermal  noise  process  from  each  element  into  quadrature  (orthogonal) 
components,  multiplying  the  quadrature  components  by  real -valued 
weighting  coefficients,  and  summing  the  weighted  products  to  form  the 
array  output.  These  operations  will  be  modeled  mathematically  by  an 
inner  product  of  two  vectors:  an  input  vector  of  complex  envelope 
signals  associated  with  the  real  IF  input  signals,  and  a  complex  weight 
vector  representing  the  array  weighting  functions.  The  input  vector 
is  defined  as  the  column  vector 


(1)  x(t)« 


X-j(t) 

ik(t) 


where  the  k^1  component,  xi^t),  is  the  complex  envelope  (or  pre-envelope) 
of  the  total  signal  plus  thermal  noise  process,  x|<(t),  applied  to  the 
quadrature  hybrid  input  from  the  k*'1  array  element.  The  pre-envelopes 
are  defined  at  the  IF  center  frequency,  fc,  and  contain  the  amplitude 
and  phase  modulation  of  the  IF  bandpass  processes,  xjjt),  i.e.. 


(2) 


xk(t)  =  Re  jxk(t)  • 
■  *  |xk(t)|- 


l 


cos[o>ct  +  Arg  xk(t)]. 


k=l  ,2,*  •  •  ,m. 


*A  common  local  oscillator  signal  is  applied  to  the  mixer  In  each 
element. 


Ist  ELEMENT 


SIGNAL 

Fig.  1.  Adaptive  array  structure. 


A  decomposition  of  the  input  vector,  Eq.  (1),  into  direct  signal, 
multipath  signal,  and  thermal  noise  components  will  be  discussed  later 
in  Section  II-B.  The  array  weight  vector  is  defined  as  a  column  vector 
of  complex -valued  weighting  coefficients. 


w^t) 

• 

^11  (t )  3  W21  ^ ) 

• 

(3)  w(t)  = 

wk(t) 

= 

• 

wlk(t)  +  j  w2k(t) 

• 

_wm(t)J 

• 

_><*>  *  i  w2m(t^_ 

where  wik(t)  and  W2k(t)  are  the  in-Dhase  and  quadrature  weights, 
respectively,  associated  with  the  k*h  array  element.  As  illustrated 
in  Fig.  2,  the  weighted  output  process  from  the  kt(l  array  element  has 
a  pre-envelope  given  by  the  product 

(4) *  Yk(t)  *  xk(t)  wlk(t)  +  [-j  xk(t)]  w2k(t) 

*  [wlk(t)  -  J  w2k (t)Jxk(t)  *  wk(t)+  xk(t). 

The  pre-envelope  of  the  array  output  process  is  equal  to  the  sunmation 
of  suboutputs,  Eq.  (4),  over  the  array  element  index  k,  or  equivalently 
by  the  inner  product 

(5)  y(t)  =  <w(t),  x(t)>  =  w(t)+  x(t) 

-  ?  Yk(t)  =  I  wk(t)+  \(t). 
k=l  K  k=l  K  * 

The  array  Is  adaptive  as  a  result  of  weighting  coefficient 
adjustment  and  control  through  feedback  from  the  array  output.  The 
adaptive  weight-control  system  In  Fig.  1  Is  designed  to  minimize  the 
difference  (error)  between  the  array  output  and  a  reference  signal. 

The  design  Is  based  on  a  least  mean  square  error  (LMSE)  performance 
criterion[l ,5,13]  and  results  In  weighting  coefficient  responses 
governed  by  the  vector  differential  equation 


•The  symbol  t  represents  the  adjoint  (complex  conjugate  transpose) 
operator.  The  adjoint  of  a  scalar  quantity  is  simply  Its  conjugate. 
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W,k ( t) 


Fig.  2.  Weighting  coefficient  models  with 

(a)  real  signals, 

(b)  complex -envelope  signals. 


(6)  —  =  a  x(t)  e(t)+  =  a  x(t)  [r(t)  -  -r(t)]+ 

[j  *  a  [xft)1  r(t)+  -  x(t)  x(t)+  w(t)]. 

OIn  this  equation  e(t)  is  the  error  signal  pre-envelope,  r(t)  is  the 
reference  signal  pre-envelope,  and  o  is  a  positive  gain  constant  con¬ 
trolling  the  speed  or  response  (or  bandwidth)  of  the  feedback  loops. 
Figure  3  illustrates  the  vector  processing  model  defined  by  Eq.  (6). 

0 


WEIGHTING  COEFFICIENT 
MULTIPLIERS 


r  (t ) 


I _ J 

ERROR  X  INPUT  MULTIPLIERS 


Fig.  3.  Adaptive  processor  model 


For  an  m-element  array,  2m  feedback  loops  are  required.  Note  in 
Eq.  (6)  that  the  error  signal's  pre-envelope  is  conjugated  implying 
that  phase  difference  terms  of  the  error-by-signal  multiplier  outputs 
are  to  be  retained  physically.  The  double  frequency  (2fc)  components 
are  assumed  to  be  eliminated  by  low-pass  zonal  filtering  of  the  error 
multiplier  outputs  (not  shown  in  Fig.  3). 
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Equation  (6)  is  a  linear  differential  equation  in  w(t)  with 
random  coefficients,  £(t)  5f(t)+,  whose  statistics  are  dependent  on  the 
input  vector.  The  expecte?  value  of  the  weight  vector  response  to  the 
input  vector  satisfies  an  equation  given  by  the  expectation  of  Eq.  (6) 
over  all  random  variables  present: 


(7) 


dn(t) 


*]  U1L 


jt —  =  r  -  K  n(t)  -  E 
dt  -x  x  —  ' 


(x(t)  x(t)+  -  K  Kw(t)  -  n(t)}~] 


where  ojt)  =  E[w(t)]  =  mean  weight  vector 

r  =  EC_x(t )  r(t)  ]  =  input-reference  crosscorrelation 
-x  vector 

Kx  =  E[x(t)  x(t)+]  =  correlation  matrix  of  the  input  vector. 

The  latter  two  terms  on  the  left-hand  side  of  Eq.  (7)  result  from  the 
expectation  of  the  coefficient  matrix  times  the  weight  vector  processes 
in  Eq.  (6).  It  will  be  necessary  to  neglect  the  contribution  of  the 
last  term  in  Eq.  (7)  in  the  analyses  to  follow.  That  is,  the  assumption 
of  small  weight  fluctuations  about  the  mean  value,  n(t),  and  negligible 
correlation  between  weight  fluctuations  and  fluctuations  in  the  in¬ 
stantaneous  coefficient  matrix  will  be  made.  This  approximation  is 
valid  when  the  response  time  of  the  feedback  loop  is  much  larger  than 
the  reciprocal  of  the  IF  bandwidth,  i.e.,  when  loop  processing  is 
narrowband  in  comparison  to  the  modulation  bandwidth  of  the  input 
signal[9,ll].  In  general,  the  loop  gain  constant  a  can  be  adjusted  to 
satisfy  this  condition. 


The  analysis  procedure  to  be  followed  consists  of  adopting 
suitable  models  for  the  structure  of  the  input  signals,  calculating 
the  correlation  matrix  Kx  and  cross correlation  vector  rx  for  these 
models,  and  solving  the  approximate  equation 


(8) 


!  cfo(t) 
a  dt 


I*  -  Kx  l(t) 


for  the  asymptotic  value  of  the  mean  weight  vector  (as  t-**).  Suitable 
measures  of  array  performance  in  steady-state  will  then  be  analyzed  as 
a  function  of  the  parameters  assumed  in  the  models. 

B.  Received  Signal  Model 

In  this  section,  explicit  representations  for  the  input  vector 
of  pre-envelopes  and  the  associated  correlation  matrix  are  developed. 
For  analysis  purposes  an  array  of  ideal  antenna  elements  is  assumed, 
and  inter-element  time  delays  are  assumed  small  in  comparison  to  the 
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reciprocal  of  the  received  signal  bandwidth.  With  the  latter  re¬ 
striction  of  narrowband  signals  compared  to  "array  bandwidth",  the 
signals  received  from  a  specified  direction  have  identical  modulations 
in  the  different  array  elements  except  for  carrier  phase  shift.  The 
direct  (1 ine-of-sight)  signal  received  by  the  array  and  the  reflected 
(multipath)  signal  received  by  the  array  are  represented  as  separate 
vector  components  of  the  input  vector,  Eq.  (1),  in  the  expansion 

(9)  l(t)  =  ^(t)  Vj  +  C2(t)  Vg  +  n(t) 

where  £,(t)  =  pre-envelope  of  direct  signal  at  the  phase  center 
of  the  array, 

£2(t)  =  pre -envelope  of  multipath  signal  at  the  phase 
center  of  the  array, 

v.  =  vector  of  rf  phase  shifts  associated  with  the  itfl 
signal  (i =1  direct;  i=2  multipath) 


design  value  of  rf  carrier  frequency 

incremental  path  delay  from  the  array  phase  center 
to  the  kth  array  element  in  the  direction -of -arrival 
of  the  ith  signal , 

vector  of  pre-envelopes,  ?k(t),  k=l,2,***,m  of 
element  thermal  noise  processes. 

The  incremental  time  delays,  T-jk*  in  the  phase-shift  vectors  Vi  are 
dependent  on  array  geometry  ana  signal  angle-of-arrival .  For  the 
special  case  of  an  m-element,  equispaced  linear  array  of  element 
separation  p,  and  a  signal  incident  at  8^  degrees  from  the  array  axis, 
the  vector  vi  is  given  by 


“rf  = 
Tik  = 

n(t)  * 


e 


(10) 


*t*< 


3  ♦  < 


**  = 


J  0 


e"J 


or 


.2  M 


.4  *i 


“Jo 


2  *i 


-J?  *i 


m-odd  _ 


m-even 


where  4>j  is  the  electrical  phase  shift  per  element  at  the  carrier 
frequency  wrf. 


<j>.  =  2n  cos  8j  ,  i — 1 ,2 . 

1  rf  1 


As  will  be  shown  later,  the  angular  separation  of  direct  and  multipath 
signals  appears  in  the  expressions  for  array  performance  as  an  inner 
product  of  the  phase-shift  vectors  V]  and  V2: 


(11) 


<», .Yj 


m 

>«  l 


^"rf^lk'W 


k=l 


For  equi spaced  linear  arrays,  the  summation  in  Eq.  (11)  is  dependent 
only  on  the  difference  in  electrical  phase  shift  per  element  (<t>-|-<t>2) , 
i  .e. , 

sin  (94 

(12)  <v1,Vo>  =  - 7, —  ;  ♦  =  [cose, -cos 6?]. 

sin(f-)  rf  1  c 

A  purely  specular  model  for  the  multipath  signal  is  adopted  in 
this  study.  The  direct  and  multipath  signal  pre-envelopes  in  Eq.  (9) 
will  be  assumed  to  be  represented  by  the  functions 

I 

I  9 

I 


03) 


^(t) 


lz(t) 


j[p(t)  +  Aw.t  +  0-.] 

5  *  ' 

j[p(t-em)  +  A^t  +  62] 


where  v's  =  rms  amplitude  of  direct  component 

*/IT  =  rms  amplitude  of  reflected  component 

p(t)  =  angle  modulation  of  direct  component 

e  =  differential  path  delay  of  reflected  component 
relative  to  direct  component 

2  =  carrier  frequency  offsets  from  the  design  value 

0,  ,  =  carrier  phase  angles  at  the  array  phase  center. 


Frequency  offsets  representing  carrier  frequency  uncertainties  or  drift 
and/or  doppler  frequency  shift  are  included  in  the  above  representations 
for  the  sake  of  completeness.  The  effects  of  these  offsets  on  array 
response  will  not  be  determined  in  the  present  study,  however,  due  to 
difficulties  encountered  in  the  analysis;  i.e.,  these  offsets  will  be 
set  equal  to  zero  later.*  All  time-independent  parameters  of  the 
signals  in  Eq.  (13)  are  assumed  to  be  non-random.  The  angle  modulation 
p(t)  will  be  modeled  as  a  random  process  which  includes  the  transmitted 
information  modulation  (data)  plus  any  subcarrier  phase  modulation 
(code)  added  for  bandspreading  purposes.  The  pre-envelope  of  the  trans¬ 
mitted  signal,  exp[j  p(t)],  will  be  defined  to  have  a  stationary  ensemble 
autocorrelation  function 

(14)  Rp(t )  =  E[eJ’  p(t)  e"J  p(t"x)]. 

A  particular  type  of  phase  modulation  of  special  interest  is  that 
employed  in  coded,  binary  communications  systems.  The  communications 
signal  is  generated  by  bi -phase  modulating  a  carrier  signal  with  the 
modulo-two  sum  of  a  random  binary  data  sequence  and  a  pseudonoise  (PN) 
code.  This  signal  will  be  modeled  as  the  product  of  two  independent 


*A  closed-form  asymptotic  solution  to  the  differential  equation,  Eq.  (8), 
has  not  been  found  for  the  case  of  non-zero  frequency  offsets. 


-  -■*+ 


random  binary  processes  with  bit  periods  of  and  A  seconds  which 
correspond,  respectively,  to  the  data  and  code  bit  durations.  The 
autocorrelation  function  of  this  signal  is  a  product  of  triangular 
functions : 


V'>  - 


(1  -  -Lf1)!1  - 111-) 


T  <  * 


T  >  A 


<4  «  V- 

The  autocorrelation  and  crosscorrelation  functions  of  the  signals 
present  at  the  receiving  array  are  calculated  from  Eq.  (13)  as  follows: 


E[^(t)  ^(t-x)+]  =  S  eJ  Rp (t ) 

...  jA(D9T 

E[?2(t)  €2(t-T)+]  =  M  e  1  Rp(t) 


E[^(t)  ?2(t-x)+]  =  m  e 


Ot -^u>2 )t  +  AwgT  +  ©i  “e2-^ 


RP(t  +  £m> 


The  crosscorrelation  function  in  the  last  expression  is  non -stationary 
when  the  differential  doppler  shift  (A«i  -  A<*>2 )  between  direct  and 
multipath  signals  is  non-zero. 

The  thermal  noise  processes  in  Eq.  (9)  are  assumed  to  Gaussian, 
stationary,  and  independent  of  each  other  and  the  received  signals. 
These  noise  processes  are  assumed  to  have  a  constant  power  spectral 
density,  N0,  over  the  noise  bandwidth,  Bif,  of  the  IF  amplifiers.  The 
correlation  matrix  associated  with  the  input  thermal  noise  vector  is 
given  by  the  diagonal  matrix 

(17)  E[n(t)  n(t)+]  =  o2  I, 


where 


0  =  Bif  No 


is  the  mean-square  amplitude  of  each  noise  process.  (I  is  the  identity 
matrix. ) 
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Consider  now  the  calculation  of  the  correlation  matrix,  Kx.  For 
narrowband  input  vectors,  as  in  Eq.  (9),  this  matrix  is  given  by 

(18)  Kx  =  E[x(t)  x(t)+] 

=  E[|?f1(t)|2]  v1v1++E[l€2(t)|2]  v^ 

+  E[C-,  (t)  £2(t)+]  v,v^  +  E[if2(t )  C1(t)+]  VgV.|  + 

+  E[n(t)  n(t)+]. 

In  this  result,  a  dyadic  matrix  v-j v-j +  is  present  corresponding  to  each 
input  signal,  and  a  pair  of  adjoint  matrices,  viVj  +  and  VjVj+,  is  present 
due  to  the  crosscorrelation  of  the  itfl  and  j™  input  signals.  For  the 
particular  model  represented  by  Eqs.  (13) -(17),  the  correlation  matrix 
has  the  form 


(19)  Kx(t)  =  S  v-jV-j1,  +  M  VgV^ 


+  «hl(o  + 'K>  W 

P  m  H-n: 


+  c2  I 


^1 


where  “12  5  Awl  " 

012  5  ®1  '  02* 


This  correlation  matrix  has  elements  which  vary  periodically  with 
time;  the  frequency  of  repetition  is  equal  to  the  differential  doppler 
frequency  shift  between  direct  and  multipath  signals. 

C.  Array  Reference  Signal  Model 

Ideally,  the  reference  signal  for  the  array  is  a  replica  of  the 
desired  signal  component  in  the  array  output.  The  error  signal  in  this 
case  contains  only  undesired  components,  and  the  feedback  processing 
acts  to  remove  these  components  from  the  array  output.  In  practice. 
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however,  an  ideal  reference  signal  cannot  be  generated  since  received 
signal  parameters  may  not  be  known  precisely  and  unknown  (data)  modulation 
may  be  present  in  the  received  signal.  Practical  techniques  for 
obtaining  the  reference  signal  from  the  array  output  signal  through 
appropriate  waveform  processing  operations  have  recently  been  developed 
for  the  case  in  which  the  desired  signal  is  a  coded,  biphase-modulated 
communication  signal [12 ,13,14].  Similar  techniques  which  use  a 
"bootstrap"  reference  loop  to  generate  an  estimate  of  the  desired 
received  signal  have  been  investigated  experimentally  for  the  case  of 
FSK  modulation [15].  Although  the  reference  signal  obtained  using  these 
approaches  is  imperfect,  due  to  practical  limitations,  there  is  con¬ 
siderable  analytical  and  experimental  evidence  which  indicates  that 
the  effect  of  these  imperfections  may  be  minimized  through  appropriate 
design[13]. 

In  the  present  study,  the  reference  signal  will  be  modeled  as 
a  non-processed  signal  which  is  independent  of  the  array  weighting 
coefficients  and  which  contains  the  angle  modulation  of  the  input 
signal  to  within  a  timing  error,  er,  and  a  carrier  frequency  offset, 

A«r: 


(20) 


_  j[p(t-e  )  +  Ao)  t  +  0  ]. 
r(t)  =  /Re  r  r  r 


The  timing  error  parameter  in  this  representation  may  be  varied  for 
perfect  time-base  alignment  with  the  direct  signal  (er=0),  or  multi  - 
path  signal  (er=em),  or  for  partial  alignment  near  these  values  of  delay. 
The  input-reference  crosscorrelation  vector  may  be  calculated  from 
Eqs.  (9),  (13),  (14),  and  (20)  as 

(21)  rx(t)  =  E[$(t)  r(t)+]  = 

. _  j[ (Ao), -Am  )t  +  (e,-e  )] 

-^Se  1  r  1  r  R  (er)  v., 


j[ (Aou-Aw  )t  +  (e?-e  )] 

♦  •«  •  Vv£m>  h 


The  components  of  this  vector  along  the  vectors  vi  and  v?  vary  peri¬ 
odically  in  time  at  unequal  rates  in  general.  Only  in  the  special 
cases  where  A»r=A»i  or  A«r=Au»  is  there  a  single  frequency  of 
periodicity  which  matches  that  of  the  correlation  matrix.  These 
choices  simplify  the  task  of  finding  the  forced  response  of  the 
system,  Eq.  (8);  however,  the  basic  problem  of  finding  solutions  to 
the  homogeneous  system  (rx=0)  containing  the  periodically  stationary 
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matrix,  Eq.  (19),  is  not  affected.  In  the  next  chapter,  closed-form 
solutions  are  presented  only  for  the  limiting  case  of  zero  frequency 
offsets : 


Aw,  =  Aa>0  =  Au»  =  0. 
i  c.  r 


III.  PERFORMANCE  ANALYSES 

A.  Mean  Weight  Vector  Response 

For  the  case  where  all  carrier  frequency  offsets  in  the  model 
are  assumed  to  equal  zero,  the  mean  weight  vector  equation,  Eq.  (8), 
assumes  the  form 


(22) 


do(t) 


1  urn 
a  dF 


-  A  v,  +  B  Vg 


“Is  W  +  M  W  +  aZI 


+  C  ^  +  C+ 


n_(t ) 


where 


A  =  SR/S  R  (e  )  e 
P  < 


J(0l"er) 


B  =  yf/R  Rp(er-em)  e 


i(e2-er) 


j012 


C  =  vSVff  Rp(€m)  e 


The  asymptotic  solution  to  Eq.  (22)  is 
dn  (t )  I 


dt 


=  0 


D.(^) 


I* 

f  *  +  2  ^ 

S  *,*,  ♦  *  *2*2  +  «  1  I 


(A  y^+B  Vg), 
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The  inverse  matrix  can  be  calculated  analytically  using  the  Woodbury 
identity  as  discussed  in  Reference  [13].  A  simpler  approach  yielding 
identical  results  is  to  substitute  into  Eq.  (22)  the  trial  solution 

dn  (t ) 

(23)  n(t)|t>-  =  a  Vq  +  b  Vg  ,  jt  =  0  . 

This  substitution  yields  a  linear  combination  of  the  vectors  vq  and 
V£  on  the  right-hand  side  of  the  equation  and  the  zero  vector  on  the 
left-hand  side.  The  vectors  V]  and  V2  are  linearly  independent  in 
general.  Algebriac  dependence  of  vq  and  V2  occurs  only  when  the  dif¬ 
ferential  delays,  t^,  in  Eq.  (11)  satisfy  a  very  special  set  of 
relationships.  Thus,  the  coefficients  multiplying  vi  and  V2  in  the 
left-hand  side  of  Eq.  (22)  may  be  set  to  zero  individually  yielding 
a  pair  of  simultaneous  algebraic  equations  in  the  unknown  coefficients 
of  Eq.  (23): 

(mS+a2+C<V2  ,Vq>)  a  +  (mC+S<v-j  ,v^> )  b  =  A 


The  effect  of  multipath  interference  is  an  alteration  of  the  relative 
magnitudes  of  weight  vector  components  along  the  vectors  V]  and  V2. 

This  is  reflected  by  the  presence  of  the  coefficients  B  and  C  in  the 
expressions  of  Eq.  (24).  These  coefficients  are  equal  to  zero  when 
the  input  interfering  signal  is  uncorrelated  with  the  desired  (direct) 
input  signal. 

B.  Output  Signal  Mean  Phase 

In  this  section  the  mean  phase  of  the  array  output  signal  is 
calculated  for  the  case  where  the  mean  weight  vector  approaches  its 
steady-state  value  in  Eq.  (23).  The  first  question  of  interest  is 
whether  or  not  the  mean  phase  approaches  the  mean  phase  of  the  reference 
signal  when  the  output  signal  consists  of  two  correlated  components. 

The  asymptotic  form  of  the  output  signal  is  given  by 


Ys(t)  =  n(»)  +  {^(t)  v,  +  C2(t)  v^} 
t** 


,  .+  j  _  j[p(t)+e1]  j[p(t-e  )+02] 

=  nWM^e  V|  +  /lie  m  L 

j[p(t>e  )+e  ] 

*  e  r  r  Q(t) 


where  Q(t) 


n(*)+v.1  r/s  e 


j(8l-0r)  f  j[p(t)"P(t_er^ 


The  two  terms  comprising  Q (t )  are  complex  constants  when  the  time 
delays  em  and  er  are  zero,  and  are  complex  random  processes  otherwise. 
The  expectation  of  Q(t)  represents  the  average  amplitude  and  phase 
(relative  to  e_)  of  that  portion  of  the  array  output  signal  which  is 
correlated  witn  the  reference  signal.  This  expectation  is  given  by 


.  -j[p(t-6  )+0  ]  , 

e|ys(0  e  r  r  -  E[Q(t)]| 

) 

■  IlMVj  Ss  e  1  Rp(er) 

+  l(-)+Vg  ^  e  M  Vem} 


VVem> 
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where  the  inner  products  involving  the  mean  weight  vector  are 
evaluated  as  shown  below. 


(29)  n(«)+v1  =  ma+  +  <v2,v-|>b+ 


=  (m2.  |<v  ,  >|2)  a**-*V) 

|D|  jD 


+  <v2»v-]>o 


2  _B_ 

|Dt 


(30) 


n(-)V-  ta2  -  l<v,,v2>|2)iM!-_A!cl 

t  *  £  Ini  ml 


+  <V-|  .V* 


2  _k_ 

!o| 


Each  inner  product  consists  of  the  sum  of  three  complex  numbers  whose 
arguments  are  given  by  one  of  two  possible  values: 


/a!  =/bV  =  er-e1 

(31)  t  /  + 

/b_  =/A_x  =  e-e2 


This  result  follows  from  the  definitions  for  the  parameters  A,  B,  and 
C  in  Eq.  (22)  when  the  complex  autocorrelation  function  Rp(t)  is  a 
real -valued  function  of  t.  This  condi ton  will  be  satisfied  for 
balanced  types  of  digital  phase  modulation  p(t)  whose  allowed  phase 
states  are  symmetrically  distributed  over  (0,  2tt)  and  occur  with 
equal  probability.  The  inner  product  of  vectors  V]  and  V2  in  Eqs. 

(29)  and  (30)  is  real-valued,  as  noted  in  Eq.  (12),  when  the  reference 
point  for  signal  representation  is  chosen  at  the  array  phase  center. 
Thus,  the  first  term  of  Eq.  (28),  the  direct  signal  contribution,  is 
comprised  of  three  components  of  which  the  first  two  are  real -valued. 
The  remaining  component  has  its  argument  along  (61-62)  and  is  given  by 


<v2,v1>a 


2  /RVfiVs  R  (  )  R  (  )  j(0Te2) 

ini  vvv  VV  e 
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The  same  conclusions  are  true  for  the  second  term  in  Eq.  (28)  —  the 
multipath  contribution  —  with  the  exception  that  the  phase  argument 
of  the  third  component  is  along  (02-0-j): 


<v1 ,v2>o 


2  R  (e  )  R  (e  — e  )  e^2’6’’ 

p  r'  p  r  m' 


|D| 


Upon  comparing  the  above  expressions,  it  can  be  concluded  that  the 
third  components  are  a  conjugate  pair;  hence,  their  sum  is  real-valued 
which  implies  that  Eq.  (28)  is  real-valued.  Figure  4  shows  a  phasor 


COMPOSITE 

OUTPUT 


Fig.  4.  Mean  amplitude  and  phase  of  direct  and 
multipath  components  in  array  output. 


diagram  of  the  components  of  Eq.  (28)  where  each  component  has  been 
rotated  counterclockwise  by  e_  degrees.  These  results  indicate  that 
the  steady-state  composite  output  signal  is  phase-aligned  with  the 
reference  signal  on  the  average  whereas  its  components,  the  direct  and 
multipath  output  signals,  are  not. 

The  composite  output  signal  in  Eq.  (26)  may  also  be  separated 
into  a  component  which  is  completely  correlated  with  the  reference 
signal  plus  a  component  which  is  uncorrelated  with  the  reference  signal 


(32) 


Y>) 


t+~ 


j[p(t-e  )  +  0  ] 

-  e  r  r  E[Q(t)] 


j[p(t-e  )  +  6  ] 


+  e 


{Q ( t )  -  E[Q(t)]>  . 


The  second  term  in  this  expression  represents  the  uncorrelated  output 
component.  The  correlated  output  component  contains  the  correlated 
"portion"  of  both  the  direct  and  multi  oath  signals  at  the  array  output 
while  the  uncorrelated  component  contains  the  remainder  of  these 
signals.  The  importance  of  this  decomposition  is  that  the  uncorrelated 
output  component  is  not  usable  for  data  detection  when  the  time -base 
of  a  demodulator  following  the  array  is  synchronized  with  the  time-base 
of  the  reference  signal.  A  discussion  of  post-array  data  demodulation 
is  given  later  in  Section  III-D. 


C.  Output  Multi  path- to-Oirect  Signal  Ratio 

One  indicator  of  array  performance  is  the  ratio  of  powers  in  the 
direct  and  multipath  signal  components  at  the  array  output.  With  the 
definition  of  power  in  a  signal  component  as  the  mean-square  amplitude 
of  its  associated  pre-envelope,  the  output  multi path-to-direct  signal 
power  ratio  is  given  by 


V  „  ElwU^UK,!2  .  „  lilt)*./ 
ps t  ElwUi^u)*,!2  s  WtiV,!2 


The  approximation  made  on  the  right-hand  side  of  this  equation  neglects 
the  effects  of  jitter  in  the  weight  vector  about  its  expected  value, 
n(t).  From  the  results  in  Eqs.  (29)  and  (30),  the  asymptotic  value 
of  this  ratio  as  the  mean  weight  vector  settles  to  its  final  value 
is 
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(34) 


PM  Ml<vl  »v2>a^A++mff^B++(rn^_  l<vi  »v2>  j2)  (SB+-A+C)  |2 
PS  S|ma2A++<v2,v1>o2B++(m2-|<v1 ,v2> |2) (MA+-B+Ct)j2 


<V,  » Vp>  nt  f  <V.  ,V„> 

1  2  ♦  Vh  --V- 

A  V  m 


m 


i-2  ,c  B'  r » 
ma  (S  — r-  -  C) 

At 


1  + 


<V2  >V1 > 
m 


!VV  U-2  (M  _  B^c+} 
nr  J  A 


The  computer  program  in  Appendix  A  was  written  to  evaluate  this 
ratio  under  the  assumption  that  the  signal  autocorrelation  was  given 
by  Eq.  (15)  with  Tb  »  A,  i.e.. 


(35) 


RpW  - 


t  ° 


It  I  <  A 


T  |  ^  A 


A  four-element  linear  array  was  assumed  in  the  program,  using  Eq.  (12) 
with  m=4.  The  results  are  shown  in  Figs.  5-11. 


Figure  5  illustrates  performance  as  a  function  of  the  multipath 
delay,  e™,  expressed  in  PN  chip  durations  (A  =  one  code  chip),  for  several 
values  of  the  angular  separation  parameter  \i>.  The  input  signals  are 
equal  in  amplitude  and  the  reference  signal  is  perfectly  synchronized 
with  the  direct  component  (er=0).  The  carrier  phase  difference  (61-62) 
is  zero  in  this  set  of  curves.  As  the  multipath  delay  increases  from 
a  small  value,  the  multipath  component  is  increasingly  rejected  by  the 
array.  For  multipath  delays  exceeding  one  PN  chip  duration,  the 
array  treats  the  multipath  signal  the  same  as  it  would  any  uncorrelated 
interference.  The  abrupt  change  in  the  slope  of  the  curves  at  the  trans¬ 
ition  point  (one  chip  delay)  is  attributed  to  the  abrupt  change  in  slope 
of  the  triangular  autocorrelation  function  at  its  base.  As  might  be 
expected,  a  multipath  signal  more  widely  separated  in  angle  from  the 
direct  signal  experiences  greater  suppression  at  the  array  output.* 


*An  angular  separation  of  signals  of  one-half  the  natural  half-power 
beamwldth  of  the  array  corresponds  to  a  value  of  of  approximately 
41  electrical  degrees. 
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ARRAY  OUTPUT  MULT  1  PATH- TO -OEStREO  SIGN 


Fig.  5.  Multipath-to-desired  signal  ratio  at  array  output. 


pm/ps  (dB 


Fig.  6.  Multipath-to-desired  signal  ratio  at  array  output. 


22 


(OB) 


Fig.  8.  Multi  path -to ^desired  signal  ratio  at  array  output. 
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Figure  6  illustrates  the  variation  in  Eq.  (34)  versus  the  carrier 
phase  difference  of  the  input  signals.  Several  curves  are  shown  for 
different  values  of  multipath  delay  and  reference  timing  error.  Note 
that  each  of  these  curves  has  a  maximum  and  a  minimum  at  the  endpoints, 
i.e.,  when  the  carrier  phase  difference  is  either  0°  or  180°.  Thus, the 
output  ratio  is  bounded  by  its  values  at  these  two  values  of  carrier 
phase  for  specified  values  of  the  other  input  parameters.  It  can 
be  shown  analytically  that  the  derivative  of  Eq.  (34)  with  respect  to 
carrier  phase  difference  is  zero  at  these  values.  The  proof  is  lengthy, 
however,  and  will  not  be  presented  here.  Although  the  results  in 
Fig.  6  apply  to  the  steady-state  behavior  of  the  output  ratio  versus 
carrier  phase  difference,  they  are  also  expected  to  be  valid  when  the 
carrier  phase  difference  (e-j -02)  is  varying  slowly  as  a  function  of 
time.  That  is,  for  very  small  differential  doppler  shift  between 
input  components  compared  to  the  reciprocal  time  constants  of  the 
feedback  loops,  the  variation  in  the  output  ratio  is  expected  to  be 
that  given  in  Fig.  6  where  the  abscissa  represents  a  time  interval 
equal  to  one-half  the  period  of  the  differential  doppler  frequency.  In 
essence,  the  feedback  loops  "track"  the  carrier  phase  increments  and 
maintain  optimum  "steady-state"  performance  at  each  instant  of  time  in 
this  case. 

The  results  in  Figs.  7  and  8  indicate  that  the  output  ratio 
is  very  sensitive  to  small  changes  in  the  timing  error  parameters.  The 
ratio  is  plotted  versus  multi path -reference  relative  timing  error;  an 
abscissa  value  of  zero  corresponds  to  perfect  synchronization  of  these 
two  signals.  The  parameter  distinguishing  the  different  curves  is  the 
timing  error  between  direct  and  reference  signals.  The  direct  component 
is  suppressed  at  the  array  output  when  its  alignment  with  the  reference 
is  relatively  poor,  e.g.,  when 


in  Fig.  7.  The  sensitivity  of  the  output  ratio  to  timing  error 
parameters  is  more  extreme  when  the  carrier  phase  difference  is  180° 
as  indicated  by  the  results  in  Fig.  8.  In  Figs.  9  and  10,  the  ampli¬ 
tude  of  the  multipath  component  has  been  reduced  by  3  dB  at  the  array 
input.  The  curves  in  these  figures  have  the  same  basic  shape  as 
those  in  Figs.  7  and  8  but  are  more  closely  spaced.  The  same  type  of 
behavior  is  noted  in  Figs.  10  and  11  where  input  signal  powers  have 
been  reduced  by  10  dB  relative  to  thermal  noise  power. 


Usable  Output  Signal-to-Noise  Ratio 


The  results  of  the  preceding  section  show  that  the  multipath- 
to-direct  signal  ratio  at  the  output  of  the  array  varies  considerably 
as  a  function  of  multipath  delay  and  reference  synchronization  error. 


A  very  large  value  of  this  ratio  --  occurring  when  the  multipath  com¬ 
ponent  is  more  closely  synchronized  with  the  reference  than  the  direct 
component  —  is  not  necessarily  indicative  of  poor  performance,  however, 
since  a  portion  of  the  specular  multipath  signal  at  the  array  output 
is  usable  for  data-detection  purposes.  The  adaptive  array  acts  as  a 
diversity  receiver  adjusting  the  gains  of  the  two  output  components 
according  to  their  input  amplitudes  and  relative  synchronization  with 
the  reference  signal.  A  meaningful  array  performance  measure,  therefore, 
should  reflect  the  characteristics  of  this  diversity  by  indicating  the 
proportion  of  total  power  at  the  array  output  which  is  useful  in  a 
post-array  demodulator.  To  this  end,  the  ratio  of  usable  signal 
power  to  noise  power  at  the  array  output  in  a  bandwidth  equal  to  the 
spectral  width  of  the  data  modulation  is  determined  in  this  section. 
Usable  signal  power  is  defined  here  as  the  power  in  the  signal  component 
which  is  correlated  with  the  reference  signal.  Noise  power  includes 
the  power  in  the  uncorrelated  signal  component  and  the  power  in  the 
thermal  noise  components  which  falls  within  the  data  bandwidth. 

As  a  first  step  in  the  calculation,  we  consider  the  hypothetical 
processing  operation  on  the  array  output  signal  indicated  in  Fig.  12. 


Fig.  12.  Complex-envelope  model  of  a 
post-array  demodulator. 


The  array  output  is  multiplied  by  the  reference  signal  and  the  product 
is  averaged  in  a  low-pass  filter  whose  bandwidth  is  matched  to  the 
data  rate.  This  operation  is  very  similar  analytically  to  coherent 
detection  of  a  coded  communication  signal.  The  essential  difference 
is  that  the  multiplier  reference  signal  in  Fig.  12  contains  both  code 
and  data  modulations  whereas  the  derived  reference  in  a  coherent 
detector  contains  only  the  code  modulation.  Thus,  the  filter  output 


29 


in  Fig.  12  is  a  d.c.  level  (ideally,  in  the  absence  of  thermal  noise, 
timing  errors,  and  carrier  frequency  offsets)  whereas  the  output  of  a 
coherent  detector  is  the  data  modulation.  Note  that  the  thermal  noise 
component  at  the  array  output  and  also  at  the  multiplier  output 
occupies  the  full  spread  bandwidth,  Bjf,  of  the  received  signal  (see 
Eq.  (17)).  Hence,  upon  filtering  to  the  data  bandwidth,  Bd,  the  power 
in  this  component  is  reduced  by  approximately  the  bandspreading  ratio: 
Bif/Bd-  The  bandwidth  of  the  uncorrelated  signal  component  at  the 
multiplier  output  is  also  approximately  equal  to  B-jf.  It  will  be 
assumed  that  data  filtering  attenuates  the  power  in  this  signal  by 
the  bandspreading  ratio.  With  this  approximation  the  signal -to-noise 
power  ratio  at  the  data  filter  output  may  be  written  as 


(36)  SNR  =  SNR  =  - §^(t)] _ 

/*d  \  y  /  B .  \  ? 

W  (l^J(E[y2(t)]-E2[y(t)]) 

where 

(37)  y(t)  =  multiplier  output  process 

=  Re  z(t) 


(38)  z(t)  =  y(t )  l(t)+ 


=  w(t)+|^(t)v1  +  ^(tjvg  +  n(t)| 


-j[p (t-e  )+0  ] 
e 


z(t)  is  the  complex-valued,  baseband  process  at  the  multiplier  output 
associated  with  the  complex  envelope  model  of  Fig.  12.  Its  real  part, 
denoted  above  by  y(t),  is  the  real  signal  at  the  multiplier  output 
having  physical  significance: 


(39)  y(t)  -  Re  z(t)  =  |?(t)||»(t)|  cos[/r(t)  -  /*(t)]. 
The  asymptotic  mean  value  of  z(t)  is  given  by 


(40) 
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This  result  follows  from  Eqs.  (38)  and  (28)  upon  noting  that  the  thermal 
noise  component  in  the  multiplier  output  has  zero  mean  value.  Since 
the  asymptotic  mean  value  of  z(t)  was  shown  to  be  real -valued  in  the 
analysis  following  Eq.  (28),  it  follows  that  the  asymptotic  mean  value 
of  y(t)  is  also  given  by  Eq.  (28),  i.e., 

(41)  E[y(t)]|t+i>  =  Re  E[z(t)]|t^.  =  E^t)])^ 

where  E[z(")]  =  a  real  number. 

Therefore,  the  square  of  Eq.  (28)  gives  the  usable  signal  power  in  y(t): 


(42) 


?[y(t)]|t^,  =  S |n.(-)+v1 12  Rp2Ur)  +  M|n(")+v2|2  Rp2(er-em) 

j(e1-e2) 


+  2/sFf  Rp(er)  Rp(£r-Em)  ReTe 


*  (n(“)+v1  )(l(“)+V2)  + 


The  mean  square  value  of  y(t),  i.e., 

(43)  E[y2(t)]  =  E[|z(t)|2  -  {Im  z(t)}2], 

does  not  appear  to  be  easily  evaluated.  An  upper  bound  on  the 
variance  of  y(t)  in  Eq.  (36)  may  be  obtained  by  neglecting  the  noise 
power  contribution  from  the  imaginary  part  of  z(t): 

E[|z(t)|2  -  E2[y(t)]  >  Var[y(t)]. 

This  upper  bound  is  actually  equal  to  the  variance  of  the  complex 
process  z(t),  in  the  limit  as  t-n»,  since 

(44)  VarCz(t)]  =  E[|z(t)  -  E[z(t)]|2] 

=  E[|z(t)|2]  -  |E[z(t)]|2 

and  the  mean  values  of  y(t)  and  z(t)  are  equal  and  real -valued  in 
this  limit  from  Eq.  (41).  Thus,  a  useful  lower  bound  on  the  signal - 
to-noise  ratio  at  the  multiplier  output  is  given  by  the  signal -to- 
noise  ratio  In  the  complex  signal  z(t): 


SNR  .  rngat  < 


'z  Var[z(t 


fcm  "  SNRy 


The  mean-square  magnitude  of  z(t)  is  equivalent  to  the  total  power 
present  at  the  array  output  as  seen  from  the  expression 

(46)  E[  jz(t)  |2)  =  E[|?(t)|2|*(t)|2]  =  E[|y(t)|2]. 

This  quantity  is  calculated  as  follows: 

(47)  E[|Y(t)|2]|t^  =  E[|w(t)+  £,(t)vq  +  ^(tjvg  +  n(t)  |2]|t 

=  S|n(“)+v-||2  +  M|n(»)+y2|2 


+  2/SM  Rp(E|n)  Re  <  e 


j(el-e2) 


(a.(-)1‘v1)(rL(-)+Y2)+ 


+  a2!n(-)|2  . 

The  variance  of  z(t)  follows  as  the  difference  in  Eqs.  (47)  and  (42). 
Upon  combining  the  results  in  Eqs.  (36)-(47),  a  lower  bound  is  obtained 
for  the  signal -to-noise  ratio  at  the  filter  output.  This  bound  is 


SNR  >  gJ-E-  SNR  = 
0  ”  Bd  2 


^-|a(-)t»,|2  Rp2(er)  t  ^5-Jn(-)+v2|2  «p2(tr-*,) 

B.A  +  2  Rp  WV'n,)Re(X> 


d  T|l(")+il  |2[1-Rp2(cr)]  +  !V|n.<”)%|2i:i-flp2<Vem>1 

0  r  0  r 

*  2  | h  T  tRp<en.)  -  Rp<er>Rp<V£r):l  Re(x) 


where 


j  (6  1  "®0  )  f  X  4. 

X  =  e  12  (IL(-)+v1  )(n_(«)+y2)+  . 


Note  that  the  sum  of  the  three  numerator  and  four  denominator  terms 
in  this  bound  is  equal  to  the  total  power  at  the  array  output,  Eq.  (47), 
normalized  by  the  element  noise  power,  a?.  In  the  absence  of  either 
multipath  interference  or  reference  timing  error  (M=er=0),  the  bound 
in  Eq.  (48)  reduces  to  the  expected  result  for  an  array  which  is 
co-phased  to  the  direct  signal,  i.e.. 


(49) 
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as  verified  using  Eqs.  (17),  (23),  and  (24)  where  the  mean  weight 
vector  is  along  the  direct  signal  vector  v-| : 
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The  ratio  in  Eq.  (49)  is  independent  of  the  bandspreading  ratio 
whereas  the  signal-to-thermal  noise  ratios  at  array  inputs  and  output 
decrease  as  the  spread  bandwidth,  B-jf,  is  increased. 

The  computer  program  in  Appendix  B  was  written  t.o  evaluate 
Eq.  (48)  for  the  case  of  a  four-element  linear  array.  The  signal 
autocorrelation  was  assumed  to  be  given  by  Eq.  (15)  and  the  bandspreading 
ratio  was  equated  to  the  ratio  of  code  and  data  rates,  i.e.. 


(50) 


The  signal-to-noise  ratio  at  the  data  filter  output  with  multipath 
absent,  Eq.  (49),  was  fixed  at  ten  decibels.  Figure  13  illustrates 
computed  results  for  the  same  input  power  levels  and  angular  separation 
as  in  Figs.  7  and  8.  The  two  curves  shown  for  each  value  of  reference 
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timing  error  represent  the  upper  and  lower  limits  of  data  filter  out¬ 
put  signal-to-noise  ratio  (OSNR)  ve**sus  carrier  phase  difference.  As 
noted  previously,  the  spread  between  these  curves  is  an  indication 
of  the  range  over  which  the  OSNR  would  vary  periodically  in  time  for 
very  small  differential  doppler  shift  between  input  components.  It 
is  apparent  in  Fig.  13  that  this  spread  is  largest  when  both  the  input 
components  are  closely  synchronized  with  the  reference  signal. 

Consider  the  top  pair  of  curves  in  Fig.  13  which  correspond  to 
perfect  reference  signal  timing.  In  the  absence  of  multipath,  the  OSNR 
would  equal  10  dB;  with  perfectly -timed  multipath  present,  the  OSNR 
varies  between  8.4  dB  and  15.2  dB  as  the  carrier  phase  difference  varies 
between  180°  and  0°.  When  the  multipath  signal  is  delayed  by  exactly 
one  PN  chip  relative  to  direct  and  reference  signals  (i.e.,  an  abscissa 
value  of  plus  one),  the  OSNR  is  8.2  dB  and  is  independent  of  carrier 
phase  difference.  The  loss  of  1.8  dB  relative  to  the  multipath- 
absent  case  is  the  result  of  two  effects.  One,  of  course,  is  the 
increase  in  data  filter  output  noise  resulting  from  the  presence  of  the 
uncorrelated  multipath  component.  The  other  is  due  tc  a  limitation  in 
the  spatial  resolution  capabilities  of  the  array,  i.e.,  the  array  cannot 
reject  an  uncorrelated  interfering  signal  which  is  in  close  proximity 
(i|i=45°)  to  the  desired  signal  without  a  loss  in  signal -to-thermal 
noise  ratio  at  the  array  output.  It  can  be  shown  from  the  results  in 
Reference  [13]  that  the  loss  due  to  this  effect  is  equal  to  1.3  dB  for 
the  input  power  ratios  and  angular  separation  of  Fig.  13.  In  the 
absence  of  multipath,  the  array  output  signal -to-thermal  noise  ratio 
is  equal  to  6  dB;  with  equal -ampli tude,  uncorrelated  multipath  present, 
the  direct  signal -to-thermal  noise  ratio  at  the  array  output  reduces 
to  4.7  dB.  The  uncorrelated  multipath  noise  in  the  array  output 
further  reduces  the  effective  signal-to-noise  ratio  to  4.2  dB,  and  the 
processing  gain  in  the  demodulator  of  2.5  numeric  (4.0  dB)  increases 
the  ratio  to  8.2  dB  at  the  filter  output:  the  calculated  value  in  Fig.  13.* 

Consider  now  the  lower  pair  of  curves  in  Fig.  13  corresponding  to  a 
reference  timing  error  of  two-tenths  of  a  PN  chip  relative  to  the 
direct  signal.  For  uncorrelated  multipath  (i.e.,  an  abscissa  value  of 
plus  one),  the  OSNR  is  2.1  dB  which  is  6.1  dB  less  than  the  value 
obtained  above  with  zero  reference  timing  error.  This  additional 
loss  is  partially  accounted  for  by  a  reduction  in  correlated  signal 
power  in  the  array  output  in  the  amount 


*The  non-integer  bandspreading  ratio  in  Fig.  13  would  not  be  implemented 
physically.  This  value  was  chosen  for  numerical  convenience  to  obtain 
unity  array  input  signal -to-thermal  noise  ratios  in  the  program. 
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(51) 


.542 


Rp2(er)  =  Rp2(.24)  =  C(.8)(.92)]2  = 

=  -  2.66  dB. 

The  remainder  is  accounted  for  by  the  increase  in  noise  contributed 
by  the  uncorrelated  portion  of  the  direct  signal.  That  is,  the  direct 
signal -to-thermal  noise  ratio  of  4.7  dB  (3.0  numeric)  with  zero 
reference  timing  error  is  effectively  reduced  to 

<52>  ' -685  =  - '-65  dB 

when  the  timing  error  is  two-tenths  of  a  PN  chip.  The  multipath  "noise" 
in  the  array  output  further  reduces  this  ratio  to  -  1.9  dB,  and  the 
4.0  dB  processing  gain  in  the  demodulator  improves  the  ratio  to  the 
calculated  value:  2.1  dB.  In  summary,  the  ideal  OSNR  of  10  dB  with 
multipath  absent  is  decreased  by  7.9  dB  when  uncorrelated  multipath 
is  present  of  which  1.3  dB  is  accounted  for  by  angular  separation 
effects,  2.66  dB  is  due  to  correlated  signal  loss,  and  3.94  dB  is 
due  to  the  increase  in  noise  power.  The  significant  observation  to  be 
made  here  is  that  the  uncorrelated  signal  components  at  the  array 
output  are  larger  than  the  thermal  noise  components  (compare  the 
denominator  terms  in  Eq.  (52)).  This  condition  is  primarily  the 
result  of  choosing  a  very  small  bandspreading  ratio  in  Fig.  13.  The 
bandspreading  ratio  was  increased  to  a  value  of  one-hundred  in  Fig. 

14.  With  uncorrelated  multipath  present  (abscissa  value  of  plus  one), 
the  decrease  in  OSNR  with  reference  timing  error  is  considerably 
less  than  that  in  Fig.  13.  That  is,  increasing  the  bandspreading 
ratio  tends  to  raise  the  ends  of  the  curves  in  Fig.  13.  It  is  im¬ 
portant  to  note  here  that  increasing  the  code  rate  (for  a  fixed  data 
rate)  corresponds  to  decreasing  the  code  bit  duration,  4.  Thus,  the 
improvement  cited  above  —  which  occurs  when  the  reference  timing 
error  is  maintained  a  fixed  percentage  of  the  code  bit  duration  —  is 
realized  only  when  reference  timing  error  also  decreases  on  an  absolute 
time  scale. 

The  values  at  the  centers  of  the  top  pair  of  curves  in  Fig.  14 
(where  er=0)  are  the  same  as  those  in  Fig.  13,  i.e.,  the  performance 
with  perfectly  synchronized  signals  is  independent  of  the  bandspreading 
ratio.  As  the  reference  timing  error  approaches  a  PN  chip  duration, 
the  center  values  of  the  curves  in  Fig.  14  converge  to  a  common  value  of 
OSNR:  the  same  value  as  occurs  when  the  timing  error  parameters  are 

e  e  -e., 

J1  =  0  and  -I—1  =  +  1.0. 
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Fig.  14.  Array  output  ratio  of  effective  signal  power-to-total 
noise  power  (in  the  information  bandwidth). 


This  value  is  equal  to  9.8  dB.  Note  the  reversal  of  roles  here;  the 
perfectly- timed  multipath  component  is  the  desired  component  while  the 
direct  component  becomes  uncorrelated  with  the  reference  as  reference 
timing  error  increases  to  one  PN  chip.  Also  note  that  the  center 
values  of  the  dashed  curves  in  Fig.  14  decrease  monotonically  toward 
9.8  dB  with  reference  timing  error.  In  contrast,  the  center  values  of 
the  solid  curves  first  decrease  slightly  then  increase  toward  9.8  dB 
as  reference  timing  error  is  increased.  Similar  behavior  is  noted 
for  the  center  values  of  the  curves  in  Fig.  13;  these  values  converge 
to  an  OSNR  of  8.2  dB  as  reference  timing  error  approaches  one  PN  chip. 

The  results  in  Fig.  15  show  that  performance  improves  when  the 
angular  separation  between  direct  and  multipath  components  is  larger. 
The  bandspreading  ratio  is  equal  to  ten  here  and  \p  has  been  increased 
to  sixty  degrees.  For  perfect  reference  signal  timing  (er=0),  the 
OSNR  has  a  minimum  value  of  9.1  dB  as  the  multipath  delay  and  carrier 
phase  difference  parameters  are  varied. 

The  results  in  Fig.  16  illustrate  performance  behavior  when  the 
multipath  signal  is  weaker  than  the  direct  signal.  The  parameters  are 
the  same  as  in  Fig.  13  except  for  a  3  dB  reduction  in  the  multipath 
signal.  The  performance  curves  are  nearly  identical  to  those  in 
Fig.  13  except  for  the  center  region  which  corresponds  to  the  condition 
where  the  multipath  signal  is  nearly  aligned  with  the  reference.  The 
curves  in  Fig.  16  are  depressed  in  this  region  relative  to  those  in 
Fig.  13.  This  behavior  reflects  the  lesser  utility  of  the  multipath 
signal  for  data  demodulation  purposes  when  its  power  is  decreased 
relative  to  thermal  noise  power. 

The  results  in  Figs.  13-16  represent  a  partial  description  of 
array  performance  in  a  multipath  envirnonment.  Certain  trends  in  basic 
behavior  are  apparent  in  these  results.  An  important  question  at  this 
point  is  the  following:  What  design  restrictions  must  be  imoosed  on 
bandspreading  ratio,  reference  timing  error,  and  input  angular 
separation  to  insure  that  worst-case  performance  is  acceptable  regard¬ 
less  of  the  multipath  signal  amplitude,  delay,  or  carrier  phase  angle? 
If  acceptable  performance  is  represented,  for  example,  by  a  decrease 
in  the  OSNR  ratio  of  no  more  than  one  decibel  from  its  ideal  value  of 
10  dB  for  no  more  than  fifty  percent  of  the  time,  what  restrictions 
must  be  imposed?  A  complete  answer  to  this  question  must  take  into 
account  other  effects  besides  those  studied  here;  however,  the  results 
in  Figs.  13-16  indicate  that  the  inequalities 
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Fig.  15.  Array  output  ratio  of  effective  signal  power-to-total 
noise  power  (in  the  information  bandwidth). 


Fig.  16.  Array  output  ratio  of  effective  signal  power-to-total 
noise  power  (in  the  information  bandwidth). 
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would  most  likely  need  to  be  satisfied  as  minimum  requirements.  The 
shape  of  the  performance  curves  might  also  vary  somewhat  as  a  function 
of  array  geometry  and  number  of  array  elements.  Additional  computer 
evaluations  are  needed  to  resolve  these  questions.  Additional  analyses 
are  also  needed  to  determine  the  multipath  response  characteristics  of 
delay-lock  loop  tracking  receivers [14].  These  receivers  are  used  to 
estimate  code  timing  (er)  from  the  array  output  signal  in  practical 
implementations.  For  many  applications,  the  locked  local  code  in  the 
delay-lock  receiver  is  used  in  generating  the  reference  signal  for  the 
array.  The  results  of  this  study  indicating  the  timing  error 
characteristics  of  the  locked  local  code  versus  input  parameters  should 
then  be  merged  with  results  similar  to  those  in  Figs.  13-16  to 
determine  steady-state  operating  characteristics  of  the  composite 
system.  A  final  refinement  should  also  include  an  investigation  of 
the  characteristics  of  processed-reference  signals,  i.e.,  of  "bootstrap" 
reference-generation  loops,  with  multipath  of  arbitrary  differential 
doppler  shift. 


IV.  CONCLUSIONS 

The  steady-state  response  of  an  adaptive  antenna  array  to  a 
received  signal  which  contains  a  specular  multipath  component  has  been 
examined.  It  was  found  that  closed-form  solutions  for  the  response 
could  be  obtained  when  the  difference  in  carrier  frequencies  of  the 
direct  and  reflected  components  was  small  and  the  carrier  frequency 
offset  of  the  array  reference  signal  was  negligible.  The  results 
show  that  the  array  acts  like  a  diversity  receiver  which  weights  the 
input  components  in  accordance  with  their  relative  correlation  with 
the  array  reference  signal.  This  correlation  depends  on  the  time-delay 
of  the  input  components  relative  to  the  time-base  of  the  reference 
signal  and  on  the  structure  and  banctoidth  of  modulation  contained  in 
these  signals. 

It  is  shown  that  the  array  adjusts  the  amplitudes  and  carrier 
phases  of  the  output  signal  components  in  such  a  manner  that  the 
carrier  phase  of  their  resultant  sum  is  aligned  on  the  average  with 
the  carrier  phase  of  the  array  reference  signal. 

The  results  also  show  that  not  all  of  the  power  in  the  output 
signal  components  is  usable  for  data  demodulation  purposes  as  a 
consequence  of  delay  differences  in  their  modulations.  The  ratio  of 
usable  signal  power  to  unusable  signal  power  plus  thermal  noise  power 
at  the  output  of  a  coherent  detector  Is  calculated.  The  detector 
bandwidth  Is  matched  to  the  data  rate,  and  Its  reference  signal  Is 
matched  In  delay  and  carrier  phase  to  the  array  reference  signal. 
Numerical  evaluations  of  this  performance  ratio  are  presented  for  the 
case  of  a  four-element  linear  arrey  where  the  desired  Input  signal  Is 
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a  bi-phase  modulated,  pseudonoise-coded,  communication  signal.  These 
results  show  that  ideal  performance  is  obtained  when  the  multipath 
component  is  absent  and  the  array  reference  signal  is  perfectly 
synchronized  with  the  direct  component.  Performance  degrades  when 
multipath  is  present  and  the  reference  signal  is  imperfectly  aligned 
with  either  input  component.  The  amount  of  degradation  is  dependent 
on  several  parameters.  It  is  found  that  worst-case  performance  with 
respect  to  the  relative  carrier  phases  of  the  input  components  occurs 
when  the  carrier  phase  difference  (at  the  array  phase  center)  is  equal 
to  180°.  It  is  also  shown  that  increasing  the  bandspreading  ratio, 
i.e.,  the  ratio  of  code-to-data  rates,  lessens  the  performance  degra¬ 
dation  provided  the  timing  error  of  the  direct  component  is  maintained 
at  a  very  small  fraction  of  the  pseudonoise  code  bit  duration.  The 
performance  degradation  is  also  less  severe  when  the  angular  separation 
of  input  components  is  increased.  In  general,  the  angular  separation 
must  exceed  one-half  the  (natural)  half-power  beamwidth  of  the 
(co-phased)  array  to  avoid  significant  degradation  for  arbitrary 
amplitude,  delay,  and  carrier  phase  of  the  multipath  component. 
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APPENDIX  A 
COMPUTER  PROGRAM  I 
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1  C  HP  2  USA  CJMFUTEn  BASIC  LANGUAGE  PROGRAM 

2  C  MULTI  PA  TH- TO -DIRECT  SIGNAL  RATIO 

3  C  M=N UMBER  OF  ARRAY  ELEMENTS 

A  C  SPIMPUT  SIGNAL  TO  NOISE  RATIO 

5  C  MPINPUT  MULTIPATH  TO  NOISE  RATIO 

6  C  PS  PANGULAR  SEPARATUM  <  DEGREES) 

7  C  E PREFERENCE  TIMING  ERROR 

8  C  E2=MULTI PATH  TIMING  ERROR 

0  C  PCARKI ER  PHASE  DIFFERENCE  (DEGREES) 

10  C  R« ARRAY  OUTPUT  RATIO 

11  10  PRINT  aMjSl*Ml> PS  P*3 

12  20  INPUT  M#  SUNUP 

13  30  PRINT 

1  A  AO  PHI NT 

Is  60  FOR  El =0  TO  .0  STEP  .1 

16  70  FOR  E2=  •  8  TO  -.8  STEP  -.8 

17  75  FOR  PO  TO  180  STEP  30 

18  80  LET  P2=3. 1A159»P/ 180 

19  85  LET  T2=3. 1A159*T/180 

20  90  LET  U=SIN(M*P2/2)/<M*SINCP2/2) ) 

21  100  LET  V1-1-V*V/ 

22  110  LET  QP  1 -A3S( El ) 

23  120  LET  Q2=  1  - A13S<  F.2 ) 

2A  130  LET  Q3= 1 - ABSC  F.2- El  ) 

25  1  AO  IF  Q3>1 .0000OE-03  THEN  160 

26  150  LET  Q3=0 

27  160  LET  Q=G3/G1 

28  170  LET  W*SQR(M1/S1) 

29  175  LET  NO=Q* W* ( 1 +  M*S 1  * 0 l ) -M* W*S l * Q2* 0 1 

30  180  LET  N=<U+NG*CJS<T2>  ) * < V+NO*COS( T2 >  >  ♦  ( NO*SIN  <  T2  )  >*<N0*SI.N(  T2) ) 

31  185  LET  DO*Q*  W* V 

32  186  LET  DP  1  ♦M*M1  *  v/1  *  <  1 -Q*G2  > 

33  190  LET  D»<D1+DO*COS(T2>>#*2+(DO#SIN<T2>>**2 
3A  220  LET  H*A. 3A295*L0G<M 1 *N/< SI  * D) > 

35  230  PRINT  El*£2»TiK 

36  2A0  NEXT  T 

37  250  PRINT 

38  860  NEXT  F.2 

39  2  70  PRINT 

AO  280  NEXT  El 
A 1  300  END 
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APPENDIX  B 
COMPUTER  PROGRAM  II 
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FH=(Fl-l.ti)/10.0 
WrtITE(6»*0)t.R 
20  FORMAT •///♦3MEP=»f 4»1) 

00  3o 
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